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ABSTRACT 

GJ 581 is an M dwarf host of a multiplanet system. We use long-baseline interferometric 
measurements from the CHARA Array, coupled with trigonometric parallax information, to 
directly determine its physical radius to be 0.299 ± 0.010i? Q . Literature photometry data are 
used to perform spectral energy distribution fitting in order to determine GJ 581 's effective 
surface temperature T EFF = 3498 ± 56 K and its luminosity L = 0.01205 ± 0.00024L o . From 
these measurements, we recompute the location and extent of the system's habitable zone and 
conclude that two of the planets orbiting GJ 581, planets d and g, spend all or part of their orbit 
within or just on the edge of the habitable zone. 

Subject headings: infrared: stars — planetary systems — stars: fundamental parameters (radii, temper- 
atures, luminosities) — stars: individual (GJ 581) — stars: late-type — techniques: interferometric 
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The formation, evolution, and environment of 
extrasolar planets are heavily dependent on the 
astrophysical properties of their respective par- 
ent stars. In particular, the location and ex- 
tent of the system's habitable zone (HZ) are di- 
rect functions of the parent star's size and sur- 
face temperature, both of which are frequently 
determined by stellar modelling. The advent of 
long-baseline interferometry at wavelengths in the 
near-infrared or optical range has made it possi- 
ble to directly me asure these stellar a strophysical 
properties (e.g. 



Bai nes et al.l 120071 120081 [2009, 



l2010l : Ivan Belle fc von Braunll2009l and references 
therein). 

M dwarfs are popular targets for exoplanet 
searches, in part due to their low intrinsic luminos- 
ity, which makes it possible for planets with orbital 
periods of much less than a year to be located near 
or inside the system's habitable zone. GJ 581 (= 
HIP 74995, V = 10.570 H = 6.095 if = 5.837) 
is a nearby M2-3 dwarf ( Bessell 19901 : Henry et al.l 
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1994 lHawlev et all Il996t ICutri et al.l l2003h that 
hosts six currently reported planets, all detected 
by the radial yelocity techn i que a nd successively 
announced in Bonfils et al. ([20051 Neptune-mass 
planet b), lUdrv et al.l (2007 , super -Earth mass 



2009, Earth-mass 



Vogt et all (|2010L 



planets c k, d). iMavor et al 
planet e), and most recently. 
super-Earth mass planets f & g) . 

Here we present the fundamental astrophysical 
parameters of GJ 581 primarily based on interfero- 
metric observations, and comment on the resulting 
location and extent of the system's HZ as well as 
the locations of the orbiting planets with respect 
to it. Our observations are described in Sj2j and 
the calculations of the resulting stellar astrophys- 
ical parameters in 53 We discuss the HZ of the 
GJ 581 system in 2J and summarize and conclude 
in g5] 

2. Observations 
2.1. Interferometry 

Our interferometric observations of GJ 581 
were carried out using the Georgia State Uni- 
versity Center for High A ngular Resolution As 



trono my (CHARA) Array (|ten Brummelaar et al 



2005) in the context of ongoing surveys to deter- 
mine fun damental parameters o f exoplanet host- 
ing st ars ( von Braun et al. |2010h and of low-mass 
stars ( Bovaiian et al. 2010^ We used the sin- 
gle baseline CHARA Classic beam combiner in 
ii-band (X ce ntrai — 1-67 /mi) and collected data 
on 6 nights between March and June 2010. We 
observed with CHARA's longest baselines Sl/El 
(~ 330m, 19 observations) and Wl/El (~ 313m, 
8 observations). The calibrator stars HD 136257 
(V = 7.55 J? = 6.24,6> SED = 0.233 ± 0.007 milli- 
arcsecondfQ) and HD 136713 (V = 7.99, H = 
5.83,6>sed = 0.323 ± 0.023 milli-arcseconds) were 
observed along with GJ 581 in bracketed se- 
quences. These calibrator stars were chosen to 
be near-point-like sources of similar iJ-band mag- 
nitudes as GJ 581 and located at small angular 
distances from it. 

The uniform disk and limb-darkened angular 
diameters (#ud and #ld, respectively; see Tabled]) 
are found by fitting the calibrated visibility mea- 



surements (Fig. [IJ to the respective functions for 
each relation. These formulae may be described 
as n* h -order Bessel functions that are dependent 
on the angular diameter of the star, the projected 
distance between the two telescopes and the wave- 
length of observation (see equations 2 and 4 of 
lHanburv Brown et aT1ll974l ). We use the H-band 
limb-darkening coefficient (hh — 0.358) calculated 
from the PHOENIX model atmospheres code for 
the corresponding value of Jeff = 3500 K and 
logg = 5.0, as tabulated in IClaretl (|2000TP1 . The 
solutions and their respective errors are computed 
using MPFIT, a non-linear leas t-squares fitting 
routine in IDL (jMarkwardtH2009r i. For this solu- 
tion we find xleduced = 0.65, implying that the 
standard CHARA reduction pipeline may have 
overestimated our measurement errors. Thus the 
results presented here assume a true Xreduced = 1 
to remedy effects of the overestimated uncertain- 
ties. Our results are shown in Table [T] 

Finally, our interferometric data allowed for a 
search for separated fringe packets to comment on 
the potential existence of a previously unresolved 
stellar comp anion in a near face-on orbit (see 
discu ssion in Baines et al. 2010t iFarrington et al. 
20101) . which would influence the interpretation 



of the RV curves. The fringe packets from two 
stellar sources would overlap for small projected 
separations of ^0.5-5 milli-arcseconds (mas), pro- 
ducing a periodic modulatio n in the angular di - 
ameter visibility curve (e.g. Baines et al.ll201oT) . 
The amplitude of this modulation is set by the 
flux ratio of the binary components. Based on the 
range of visibility residuals in our angular diam- 
eter fit, we can rule out companions with mag- 
nitude differences SH < 2.4 mag. The fact that 
we did not find any evidence of a secondary fringe 
packet (and thus stellar-mass companion at low 
inclination angles) reinforces the dynamical stabil- 
ity argumen t of coplanar orbits in th is system as 
descri bed in IMavor et al.l (|2009l ) and iBeust et al 



(2008) by providing an empirical lower-limit con- 
straint on the orbital inclinations of the planetary 
orbits. 



1 ^SED corresponds to the estimated angular diameter based 
on spectral energy distribution fitting. 



2 We note that an uncertainty in an assumed parameter will 
only modestly influence the value of the coefficient, and 
consequentially the value of $ld- F° r example, changing 
the model Teff by ±300 K will change the resulting di- 
ameter by only a few tenths of a percent, well within our 
error budget. 
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2.2. Astrometry 

Translating the measured angular radius of 
GJ 581 into a linear radius depends on its dis- 
tance. Trigonometric parallax measurements are 
available for GJ 581 in The Gen eral Catalogue of 
Trigo nometric Stellar Parallaxes (jvan Altena et al. 
19951 TTtrig = 157.9 ± 5.6 mas), from van 
Leeuwen's reduction o f the HIPPARCOS space 
mission observations ( van Leeuwenl l20Q8i i^trig 
= 160.91 ± 2.62 mas), and from the RECONS 
group3 of n trig = 160.01 ± 1.61 mafl The 
weighted-mean parallax of the three independent 
measurements for the GJ 581 system is 160.12 ± 
1.33 mas, driven primarily by the improved RE- 
CONS parallax, and is used here in the calculation 
of GJ 581's physical size. 

3. Fundamental Astrophysical Parameters 
of the Host Star GJ 581 

A summary of our results of the fundamental 
parameters of GJ 581 is given in Table [TJ We de- 
scribe below how individual values were obtained. 

Based on GJ 581's it trig and interferometri- 
cally measured angular diameter 6*ld, we cal- 
culate its linear radius to be R = 0.299 ± 
O.OlO-B fTv Following t h e pro cedure outlined in 



§3.1 of Ivan Belle et al.1 (|2007l ). we produce a fit 
of the stellar spectral energy distribution (SED) 
to literature UBVRIJHK photometry. We ob- 
tain a value of GJ 581's stellar bolometric flux 
F BO L = (9.92 ± 0.104) x 10" 9 erg cm " 2 s" 1 
and luminosity L = 0.01205 ± O.OOO24L . Pho- 
tometry for our SED fit from the following pa- 
pers was used al on g with spe c tral template s 



from Pickles 



Koen et all (1200 



JHK S ) 



(2007) 



g witn spe ctral templates 
Weid (Il993l ). iBessell (<1990h . 
Cutri et al.1 (120031 2MASS 



Nordstrom et all (1 2 0041) . iKilkennv et~al 



Ravner et al. ( 20091 ). and Koen et al 



20101 ). The transformation from the SAAO pho- 



tometry sy^tem_to2MASS was based on the equa- 
tions in lKoen et ail (|2007l ). We show our SED fit 
in Fig. H 
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4 This supersedes the value reported in iJao et all | |2005| ). 
where details about the RECONS astrometry program at 
the CTIO 0.9m can be found. In essence, the improved 
parallax results from a reduction of 224 total frames of the 
GJ 581 field taken over 9.93 years, whereas the previous 
result was based on 122 frames over 2.95 years. 



It should be noted that, in order to minimize 
the \ 2 °f our SED fit, we obtain a non-zero value 
for extinction (Ay — 0.174 ± 0.021), which is un- 
expected for a star at a distance of around 6 pc. 
Forcing Ay = degrades the x 2 from around 2 to 
around 6. The associated difference in the calcu- 
lated value of Fbol is around 6.7%, correspond- 
ing to a 1.6% effect in the calculation of Teff (see 
Equation [T] below) in the sense that an Ay > 
produces a slightly higher Teff- We postulate 
that the extinction could be circumstellar in ori- 
Equivalent SEP fitting using NextG en syn- 
thetic templates ( Hauschildt et all 19991 ) instead 



of the empirical Pickles! (jl998FT templates yielded 



near-identical values of Tbol and Ay for the best- 
fit templates using Teff = 3300 K, log g = 5.0, 
[Fe/H] = -0.5. 

The effective temperature Teff is calculated 
based on the re-arranged form of the Stephan- 
Boltzmann equation 



Teff(K) = 2341(F BO l/0ld)' 



(1) 



where Tbol is in units of 10 -8 erg cm~ 2 s _1 
and the angular diameter #ld is in units of mas. 
We calculate GJ 581's effective temperature to be 
Teff = 3498 ± 56 K. 

We derive a mass estimate of M ~ O.3M 
from the if-band M ass-Luminosity relation in 
Delfosse et all (|2000l) . Finally, GJ 581's mass and 
radius produce log g — 4.96 ± 0.08, assuming a 
mass uncertainty of 0.05 M & . See Table [1] for a 
summary of our results. 

The most recent values for GJ 581's astrophys- 
ical properties before the ones presented h ere are 
the ones calculated in Bo nfils et al. I (120051) . They 
were subsequently adopted by the papers docu- 
menting the disco veries of the othe r five planets 
orbit i ng GJ 581 (lUdrv et all l2007t iMavor et al ' 



20091 : IVogt et al.ll201tl) . Below we compare these 
values to ours. 

Bonfils et al.l (120051) use a l/-band b olomet- 

" (|l998l) to de- 



ric correction from iDelfosse et al 
rive L = 0.013L Q 



which constitutes a 4er in- 



5 Both lK6spal et al.1 {20091 ) and lBrvden et al.1 {20091 ) report a 
70 fj,m excess in the analysis of Spitzer MIPS archival data 
of around three times the value expected from extrapo- 
lating the stellar SE D to wavelengths red ward of 2MASS 
photometry, although Brvdcn et al. (2009) classify the sig- 
nificance of their detection as marginal. 
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crease with respect to our value, using only our 
error estimates; theirs are not given. This dis- 
crepancy decreases only marginally (to 3.6a) if 
we set equal the two slightly different values for 
TTtria- The radius of GJ 581 was calculated in 



Bonfils et al.l (]2005l ) by applying their estimated 
mass, M = 0.31M W , to the Ma ss-Radius rela- 
tion in IChabrier fc Baraffd <|2000h . We find that 
their resulting value of 0.29 Rq is ~lc (or ~3.0%) 
smaller than our dire ct measurement. Fin ally, the 
estimate of Teff in Bonfils et al. ( 20051) . based 
on their calculated stellar radius and luminosity, 
is 3622 K, which is 3.5% (2.2c) higher than our 
value. 

Interestingly, some previous estimates of GJ 581's 
?eff yield results with lower temperature esti- 
mates. The t emperat u re cal ibration for M dwarfs 
presented in iBessell (1995) implies a value for 
GJ 581's T EFF of 3320 K Furthermore, two differ- 
ent approaches used in ICasagrande et al~ (l2008h 
result in Teff=3320 K and Teff=3300 K, respec- 
tively. Conversely, the (independently calculated) 
e mpirical relation between Tef f an d (V — K)q 
in Ivan Belle fc von Braur] (|2009L equation 2) pro- 
duces an effective temperature of GJ 581 of just 
over 3500 K. 

4. The Habitable Zone in the GJ 581 Sys- 
tem 



A system's traditional HZ is defined as the 
range of circumstellar distances from a star within 
which a planet could have liquid water on its sur- 
face, given a dense enough atmosphere. The vari- 
ous criter ia for defining the H Z are described in 
detail in iKasting et al.l (Il993 ) and further gen- 
eralized by lUnderwood et al.l ([20031 ). For the 
GJ 581 system, the HZ b o undaries were calcu - 
late d bvlUdrv et al.l (120071). iMavor et al.l (|2009l ). 
and lVogt et al.l ( 20101 ) in the context of the respec- 
tive detections of additional planets in the GJ 581 
system. A detailed analysis of the habitablility 
of the s ystem under variou s assumptions can be 
found in lSelsis et all (|2007f ). 

Our recalculation of the boundaries of the 
GJ 581's HZ based on our directly determined 
ho st star properties i s bas ed on the equations 
Underwood et al. (2003) and IJones fc Sleep! 



temperature of the host star. This results in an 
inner boundary of 0.11 astronomical units (AU) 
and an outer boundary of slightly above 0.21 AU. 
To determine which of GJ 581's planets are located 
in the HZ, we calculate their equilibri um temper 



atures T eq following the methods of ISelsis et al 



(2007): 



eq ~ 



5(1 - A) 



(2) 



(|2010T ) that relate the radii of the inner and outer 
edges of the HZ to the luminosity and effective 



where S is the stellar energy flux, A is the Bond 
albedo, and a is the Stefan-Boltzmann constant. 
The redistribution factor / is determined by the 
efficiency of atmospheric heat redistribution effi- 
ciency and is set to 2 for a hot dayside and to 4 
for even heat redistribution. Table [2] shows the 
calculated equilibrium temperatures for the plan- 
ets in the GJ 581 system assuming an Earth Bond 
albedo of 0.29 in each case. We look at two differ- 
ent scenarios of the geometry of the GJ 581 system 
below. 

Scenario 1: W e take orbital parameters from 
Mayor et al] (|2009l) . Figure [3] (left panel) depicts 
this scenario of the geometry of the GJ 581 system 
with its HZ shown as a gray shaded region. Plan- 
ets d spends part of its orbit in the HZ. It should 
be noted, that due to the non-zero orbital eccen- 
tricity, T eq is a function of time (or phase angle). 
Thus, the T[ q =A for planet d varies from 229 K 
at periastron to 154 K at apastron, causing it to 
periodically dip into and out of the HZ. 

S cenario 2: We t ake all orbital parameters 
from Vogt et al. ( 2010l). motivated by the discus- 
sion in Anglad a-Escude et al. (2010), which points 
out that the RV signal of planet d's eccentric or- 
bit could, in fact, be due to two circular orbits 
of planets d and g. The geome t ry of this sce- 
nario is illustrated in Vogt et ail (|2010l . figure 6 
and table 2) and shown here in the right panel of 
Figure [3J In this scenario, planet g spends all its 
their circular orbit in the HZ with a semi-major 
axis of a g = 0.15 AU, while planet d is right on 
the outer edge of the HZ with a d = 0.22 AU. 
For planet g, an even heat redistribution results 
in T eq = 222 K w hich is sligh t ly les s than the 
value calculated by IVogt et al.l (j2010l) as antici- 
pated from our smaller value for GJ 581's luminos- 
ity. We note that T eq = 222 K is below the freez- 
ing temperature of water, but does not take into 
account the greenhouse effect heating due to an at- 
mosphere (for comparison: Earth's T eq = 255 K). 
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Sec section 6 of Vogt et al. I (l2010h and particularly 



Wordsworth et all (|2010j ) for more detailed elabo- 
rations. Planet d's T eq = 181 K for / = 4 when 
using its semi-major axis in Equation [2] 

5. Summary and Conclusion 

In this paper, we present directly determined 
astrophysical parameters of the late-type, multi- 
planet host star GJ 581, observed as part of our 
ongoing interferometric survey with the CHARA 
Array of KM dwarfs with and without planets. We 
find small but significant differences between liter- 
ature values and our empirical results for GJ 581, 
which are explicitly shown in £[3] and Table Q] 
We measure a stellar radius that is larger than 
predicted by theory-based calibrations (3%) but 
a slightly lower than predicted Teff (3.5%), re- 
sulting in our calculated luminosity's being lower 
than the currently used literature value (by ~7%). 
Thus, we confirm the results in iBovaiian et al. 
( 20101 ) that theoretically calculated radii of M 



dwarfs are smaller than the directly measured 
counterparts. 

We recalculate the boundaries of the system's 
HZ based on our new results. An analysis of the 
equilibrium temperatures of the six planets orbit- 
ing GJ 581 confirms that planet g spends all of its 
orbit inside the HZ. Depending on the assumed 
eccentricity of planet d's orbit, it periodically dips 
into and out of the HZ or spends all of its circu- 
lar orbit on its outer edge (gj). The presence of 
any moderately dense atmosphere around planets 
d and particularly g could provide sufficient green- 
house heating of the planetary surface tempera - 



tures (|Selsis et al.ll2007t IWordsworth et al1l2010h 



such that any existing water could be in liquid 
form. 

Planetary characterization is playing an in- 
creasingly important role in exoplanet research. 
With the exception of directly imaged planets, the 
large majority of astrophysical exoplanet param- 
eters are functions of stellar parameters, which 
in turn are, by necessity, calculated on the ba- 
sis of stellar models. Particularly in the low-mass 
regime, however, these stellar models lack obser- 
vational constraints, leading, not necessarily by 
fault of the models themselves, to uncertainties 
and sometimes systematic discrepancies between 
theoretical and observational results, as we present 



in this paper. It is thus hard to overstate the 
importance of "understanding the parent stars" . 
With recent and ongoing improvements in both 
sensitivity and spatial resolution of near-infrared 
and optical interferometric data quality, we are 
able to ameliorate this situation and provide firm, 
direct measurements of stellar radii and effective 
temperatures in the low-mass regime. 
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Fig. 1. — Calibrated visibility observations along 
with the limb-darkened angular diameter fit for 
GJ 581. For details, see JTQ 
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Fig. 2. — Spectral energy distribution fit for 
GJ 581. Literature photometry data are plotted 
on top of the spectral template. Horizontal error 
bars correspond to filter bandwidths. Crosses indi- 
cate flux levels of the spectral template integrated 
over the filter bandwidth. The fractional residuals 
are given in the bottom panel, along with scaled 
error bars in the respective photometry datapoint. 
The large J7-band error bar is due to uncertainties 
associated with literature photometry. For details, 
see gland Table ED 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 3. — A top-down view of the GJ 581 system. The habitable zone is indicated by th e gray shaded regio n 
with calculated boundaries of 0.11 and 0.21 AU. Left panel: scenario 1 adopted from iMavor et al.l rt2009) . 
Planet d spends part of its elliptical orbit in the HZ. Right panel: scenario 2 adopted from lVogt et al.l (|2010l ). 
Planet g spends all its orbits inside the HZ, which planet d orbits right on its outer edge. Planet f (a = 
0.758 AU) is not shown for purpose of clarity. For details, see §4] and Table [2 
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Table 1 
Stellar Properties of GJ 581 



Parameter Value Reference 



Spectral Type . . M2.5V Henry et al. (19941 

V-K 4.733 ±0.083 Bessel (19901 : Cutri et al. (20031 

[Fe/H] -0.135 Average a 

9 VD (mas) 0.433 ± 0.014 this work 

6> LD (mas) 0.446 ± 0.014 this work 

Parallax (mas) .. 160.12 ±1.33 this work 

Radius (i? ) .... 0.299 ± 0.010 this work 

Luminosity (L Q ) 0.01205 ± 0.00024 this work 

Teff (K) 3498 ± 56 this work 

logg b 4.96 ±0.08 this work 



Note. — Properties of GJ 581. For details, see $3] 

Average from iBonfils et al.l (120051); [Johnson fe Appd (|2009h ; 



Roias-Avala et al.l (|201dh ; ISchlaufman fc Laughlinl ^oibTT 



b The surface gravity logg is based on our quoted values of mass and ra- 
dius, and its uncertainty is calculated from our empirically determined radius 
uncertainty. We assume oM = O.O5M0 ~ 16%, which represents a slightly 
more conservative approach than the (aM) max = 10 % requirement in the 
calculation of the M-L relation in Delfosse et al. I (|2000h . 
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Table 2 

Equilibrium temperatures for the GJ 581 system planets. 



Planet 


a (AU) 


Tit ( K ) 


Ti= 2 ( K ) 


b 


0.041 


418 ±3 


498 ±3 


c 


0.070 


320 ±2 


381 ±2 


d 


0.220 


181 ± 1 


215 ±2 


c 


0.030 


489 ±2 


582 ±3 


f 


0.758 


97 ±1 


116 ±1 


g 


0.146 


222 ±2 


264 ±2 



Note. — Equilibrium temperature s based 



on the equations in ISelsis et all (|2007l ) 



and assuming o rbital element values 
from Mayor et all ( 20091 planets b-e) and 



Vogt et all (|20lfl planets f and g). / = 4 
implies perfect energy redistribution efficiency 
on the planetary surface, / = 2 means no 
energy redistribution. For details, see §Q and 
Fig. [3] Note that the quoted uncertainties 
are based only on the uncertainties in the 
calculation of GJ 581's luminosity, and not on 
uncertainties in the orbital elements. As such, 
they should be regarded as lower limits. 
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